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Abstract 

In  this  note  a  simple  algorithm,  which  uses  the  phase  information  from  the 
geometrical  optics  limit,  to  construct  the  shape  of  an  object  is  presented. 
Essentially,  the  phase  in  the  back  scattered  direction  determines  the  equation 
of  the  tangent  plane  at  the  unique,  but  unknown,  specular  point.  This  plane 
depends  upon  the  two  spherical  angles  8  and  4i ,  which  describe  the  incident 
wave  direction.  The  observation  that  the  tangent  plane  envelopes  the  obstacle 
as  8  and  ^  are  varied,  allows  the  derivation  of  an  explicit  formula  for  the 
equation  of  the  surface  in  terms  of  the  measured  scattered  phase.  An 
analogous  two-dimensional  formula  is  also  presented. 
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I.  Introduction 


The  field  of  inverse  scattering  is  concerned  with  the  determination  of  an 
object's  shape  from  scatterint,  data.  In  a  fundamental  paper  [1]  Keller  used 
the  theory  of  geometrical  optics  to  formulate  and  to  partially  solve  this 
problem  for  convex  objects.  Specifically,  he  showed  that  the  amplitude  of  the 
scattered  field  in  the  backscattered  direction  was  related  to  the  curvature  of 
the  obstacle  at  a  unique  specular  point  in  two  dimensions,  and  to  the  Gaussian 
curvature  in  three  dimensions.  In  the  former  case  he  was  able  to  determine  the 
obstacle's  shape  to  within  a  rigid  translation.  He  obtained  similar  results 
for  the  later  case  when  the  object  was  a  surface  of  revolution  and  the 
incident  wave  propagated  along  the  axis  of  rotation.  In  general,  he  observed 
that  the  problem  became  equivalent  to  Minkowski's  problem  whose  solution 
required  solving  a  nonlinear  partial  differential  equation.  Weiss  used 
Keller's  results  on  several  specific  examples  and  studied  the  effects  of 
noisy  data  and  of  relaxing  the  convexity  restrictions  [2].  Neither  author 
used  the  phase  information  from  the  geometrical  optics  approximation. 

Majda  and  Ralston  [3]  studied  the  direct  scattering  problem  and  found 
that  the  argument  of  the  determinant  of  the  scattering  matrix  s(A)  gave 
information  about  the  target  shape  in  the  high  frequency  limit.  Specifically 
they  derived  the  formula 

if 

s(A)  -  _  V  -  _  S  -  _  H  da  +  0(1),  A  -  ® 

32ii  8  384  aJ 

where  A  is  the  wavelength  of  the  incident  radiation,  V  is  the  volume  of  the 
target,  S  is  its  surface  area,  H  is  its  mean  curvature,  and  d  denotes  -he 
surface  of  the  target.  Using  a  result  of  Melrose  [4]  they  showed  that  s(A) 
was  related  to  the  differential  cross  section  of  the  scattered  field  in  the 
forward  scattered  direction. 


-2- 


In  this  note  we  give  a  simple  algorithm  which  uses  the  phase  information 
from  the  geometrical  optics  limit  to  construct  the  shape  of  a  convex  object. 
Essentially,  the  phase  in  the  back  scattered  direction  determines  the  equation 
of  the  tangent  plane  at  the  unique,  but  unknown,  specular  point.  This  plane 
depends  upon  the  two  spherical  angles  %  and  4>,  which  describe  the  incident 
wave  direction.  By  observing  that  the  tangent  plane  envelopes  the  obstacle  as 
0  and  4>  are  varied,  we  obtain  an  explicit  formula  for  the  equation  of  the 
surface  in  terms  of  the  measured  scattered  phase. 

II.  The  Inverse  Algorithm 

A  plane  acoustic  wave  insonifies  a  convex  target  n  which  can  either  be 
soft,  hard,  or  have  a  surface  impedance.  The  incident  wave  vector  is  given  by 

T 

Vj  =  (cos4>  sin5,sini^  sintf,  -costf)  where  ^  and  6  are  the  spherical  angles 

measured  from  a  fixed  origin  0  which  is  either  Inside  or  outside  of  fl;  see 
Figure  1.  In  the  high  frequency  limit  the  scattered  field  in  the 
backscattered  direction  is  given  by  [1,5) 

Ug  -  A  exp(-2ikXQ*  e^^ ,  as  r  -  7  x^  +  y^  +  -►  <».  (1) 

J~^  ^ 

In  this  expression  A  is  a  constant  which  depends  upon  the  particular  boundary 
condition,  Xq  is  the  position  vector  of  the  specular  point  on  the  surface  dii, 

G  is  the  Gaussian  curvature  at  this  point,  and  is  the  unit  tangent  vector 

along  the  specularly  reflected  ray.  The  location  of  the  specular  point  is 
unknown . 

We  focus  our  attention  on  the  phase  term  exp(-2ikXQ*  which  we  assume 

can  be  measured  as  a  function  of  <l>  and  0.  The  equation  of  the  tangent  plane 
at  the  specular  point  is  given  by 
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x.Vj^  -  d{<t>,e)  (2) 

where  d{<f),d)  -  Xq*  is  the  distance,  from  the  fixed  origin  0  to  the  tangent 

plane  (see  Figure  1),  which  is  known  from  the  phase  measurements.  Now 
equation  (2)  is  a  two  parameter  family  of  planes  which  envelope  the  surface 
30.  To  find  the  envelope,  i.e.  30,  we  partially  differentiate  (2)  first  with 
respect  to  ^  and  then  with  respect  to  $.  These  two  additional  equations  are 
then  appended  to  (2)  to  give  three  equations  in  three  unknowns.  They  are 

x(cos^  sin^)  +  y(sin(^  sin0)  +  z(cos6)  -  d(<^,^)  (3a) 

-x(sin.^  sin^)  +  y(cos^  sin^)  -  d^  (3b) 

x(cos^  costf)  +  y(sin^  cosO)  -  z(sintf)  -  d^  (3c) 

whore  the  subscripts  denote  partial  differentiation.  The  solution  of  this 
system  of  equations  is 

X  -  d(cos0  sin9)  -  d^(sin<^  cscff)  +  d^(cos(^  cos0)  (4a) 

y  -  d(sin^  sin?)  +  d.(cos^  cscS)  +  d.(sin^  cos^)  (4b) 

z  -  d(cos0)  -  d^(sin9)  (4c) 

which  is  just  a  parameterization  of  the  surface  30.  The  apparent  singularity 
in  (4)  at  3  -  0  is  removable,  because  d.(^,0)  -  0  as  can  be  deduced  from  (3b). 

V 

The  determination  of  the  phase  d(^,9)  is  carried  out  in  the  time  domain 
because  it  is  related  to  the  time  lag  r  between  the  Initiation  of  the 
transmitted  pulse  and  the  reception  of  the  reflected  pulse.  If  r  is  the 
distance  from  the  transmitter  to  the  origin  0,  then  an  application  of  an 
inverse  Fourier  transform  shows  the  time  between  the  leading  edges  of  the 
incident  and  reflected  pulses  to  be 

r  -  |r  +  2di<f>,9)]/c  (5) 

where  c  is  the  acoustic  wave  speed  in  the  material  surrounding  O.  For  smooth 
incident  pulses  without  a  sharp  leading  edges  this  time  is  difficult  to 
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measure  and  is  a  source  of  error.  We  shall  not  pursue  this  point  further  at 
this  time. 

III.  Examples 

In  the  first  example,  let  fl  be  a  sphere  of  radius  R  centered  at  the  point 
Xc=  (x^,y^,z^)  which  is  measured  from  the  fixed  origin  0.  The  geometrical 

optics  phase  in  the  backscattered  direction  is  -2kd  where  d  is  given  by 

=  R  +  =  R  +  x^(sin^  cos4>)  +  y^(sin(?  sin<^)  +  cosd .  (6) 

We  observe  here  that  the  phase  d  is  not  a  constant,  because  the  origin  0  is 
different  from  the  center  of  the  sphere.  Inserting  (6)  into  (4)  we  obtain 

X  =  x^  +  R(cosi^  sintf)  (7a) 

y  -  y^  +  R(sin^  sinS)  (7b) 

z  «-  +  R  cos9  .  (7c) 

which  is  just  the  equation  of  our  sphere. 

In  the  second  example,  let  0  be  the  triaxial  ellipsoid 

x^A  X  -  1  (8) 

centered  at  the  fixed  origin  0,  Here,  A  is  a  symmetric  matrix  whose 

T 

coefficients  are  to  be  determined  and  x  -  (x,y,z).  Let  P  be  the  rotation 

T 

matrix  that  diagonalizes  A,  i.e.,  PAP  D,  where  D  is  the  diagonal  matrix 
with  entries  =■  •  The  are  the  eigenvalues  of  A  and  are  the 

recipricals  of  the  axis  lengths,  and  the  columns  of  P,  ,  are  the 

corresponding  eigenvectors  which  geometrically  point  along  the  principal  axis 
of  the  ellipsoid.  In  Appendix  A  we  shown  that  the  backscattered  geometric?! 
optics  phase  is  given  by 


(9b) 


R  -  (eP)^  (eP) 

e  -  (fi.j  /  J  \  )  (9c) 

where  the  matrix  e  is  the  square  root  of  D  ^ .  According  to  the  theory 
presented  in  Section  II,  equation  (4)  gives  an  equivalent  representation  of 
(8)  using  (9).  We  do  not  demonstrate  their  equivalence  here. 

Rather,  we  pose  the  more  restricted  question:  Given  the  a  priori 
knowledge  that  fl  is  a  triaxial  ellipsoid,  how  many  measurements  must  we  make 
to  determine  the  orientation  of  the  object  and  its  principal  axis  lengths? 

In  the  low  frequency  limit  Dassios  [6j  has  shown  that  six  measurements  are  re¬ 
quired  and  we  shall  show  that  the  same  is  true  in  the  high  frequency  limit. 

Choosing  the  six  incident  wave  vectors  (1,0,0)^,  (0,1,0)^,  (0,0,1)^,  (1,1,0)^, 
T  T 

(1,0,1)  ,  and  (0,1,1)  ,  we  explicitly  determine  the  six  coefficients  of  the 
symmetric  matrix  R  from  (9a).  However,  from  the  definition  of  R,  (9b),  it  is 
T  - 1 

evident  that  PRP  =  D  ,  so  that  the  orthogonal  transformation  P  simultaneous¬ 
ly  diagonalizes  A  and  R.  Therefore,  since  R  is  known,  we  can  numerically 
compute  its  eigenvalues  and  explicitly  obtain  the  lengths  of  the  principal 
axis.  Moreover,  since  the  columns  vectors  P.  are  the  corresponding  normalized 

eigenvectors,  we  can  explicitly  determine  the  Eulerian  angles,  which  define 
the  orientation  of  n,  [7)  from  their  components. 
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IV.  Two-Dinienstonal  Result 

The  analogous  two  dimensional  result  is  derived  in  the  same  fashion  and 
we  simply  state  here  the  result 

X  -  (cos.ji)d  -  sin<jl  d'  (8a) 

y  -  (sin(^)d  +  cos<p  d'  (8b) 

where  d(4)  -  Xq* 

point,  V  •=  (cos^,sin4)  is  the  negative  of  the  incident  wave  vector,  and  the 
prime  denotes  differentiation  with  respect  to  4>. 
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specular  point,  X^,  lies  on  the 


5;i!rf.'ii  ■.  it  follox's  tnrr  (2)  onii  (,d)  that  d(4>,0)  -  1/c .  Using  the  fact  that 
P"\f'  D,  v.o  deduce  Eton  (Al>  that 

X,^  -  f'D'  '  pV,  (A2) 

I '  k 

■jtif  ■  IJ  is  ‘he  in/ei  si*  of  n.  T.drij.e  the  scalar  product  of  (A2)  with  and 

rec.  lling  the  definition  'd  d,  (2;.  v/e  obtain  the  desired  result,  (9a). 


F:)v-ro  Caption 

Fi^ute  1.  deo:netry  of  prohl'-;"  sliovirg  the  specular  point  P  ,  the  incident  and 
specviarly  r-  fleeted  •..•avo  t  limbers,  and  the  fixed  origin  0. 


